Multiscale SOC in turbulent convection 
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Abstract 

Using data obtained in a laboratory thermal convection experiment at high Rayleigh 
numbers, it is shown that the multiscaling properties of the observed mean wind 
reversals are quantitatively consistent with analogous multiscaling properties of the 
Bak-Tang-Wiesenfeld prototype model of self-organized criticality in two dimen- 



sions. 
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1 Introduction 



The existence of well-organized large-scale motion in turbulent thermal con- 
vection is an intriguing recent discovery. It has been under active investigation 
for some time, both theoretically and experimentally (see [1]-[15] and refer- 
ences therein). At least in a convection apparatus whose aspect ratio (= the 
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ratio of the diameter to the height) is unity, the large-scale motion is in the 
form of a persistent circulation (on the average), its physical extent being of 
the order of the size of the apparatus itself. This is often called the "mean 
wind" . The mean wind evolves with the Rayleigh number (which can be re- 
garded as an external tunable parameter in the flow), and seems to asymptote 
to a well-defined shape only beyond Rayleigh numbers of the order of 10^°. 
However, beginning more or less at this Rayleigh number, at least for three or 
so further decades, the mean wind undergoes abrupt and apparently irregular 
reversals of direction. These reversals are the object of the study here. 

The measurements to be analyzed were obtained at a Rayleigh number of 
1.5 X 10^^ in thermal convection occurring in a closed container of circular 
cross-section of 50 cm diameter and aspect ratio unity. The working fluid 
was cryogenic helium gas. The same data have been analyzed in the past 
[12], where the fluid dynamical origin of reversals was proposed to be the 
imbalance between buoyancy and friction (with inertia playing a secondary 
role). It was further suggested that a stochastic "change of stability" occurs 
between two metastable states corresponding to the opposite direction of the 
wind. The main physical quantity analyzed in [12] was the interval r (the 
life-time of the metastable states, or the interval between reversals). The time 
of actual switching between the two states was short on the scale of some 
average measure of r, so the wind reversals could be regarded as abrupt. The 
last property (perhaps also the precise details of reversals themselves) is most 
likely dependent on the specific boundary conditions of the experiment, but it 
is thought that the statistical properties of the duration times r are insensitive 
to these details. The analysis of [12] indicated a general dynamic mechanism 
similar to self-organized criticality (SOC) [16]-[23]. The particular emphasis 
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of this paper is the elaboration of this analogy. 

Self-organized criticality occurs through a nonlinear feedback mechanism. There 
could be many possible scenarios of SOC in our system where numerous plumes 
and jets are developed as a result of boundary layer and thermal instabilities, 
all of which are embedded in a background of strong turbulent fluctuations 
prevalent in the core. It was suggested recently for plasma turbulence (see 
Refs. [24], [25]) that instabilities governed by a threshold may lead to self- 
organized criticality by producing transport events at all scales (avalanches). 
These avalanches are due to local accumulation of energy, leading to an in- 
creasing gradient. Once the gradient exceeds an appropriate threshold, a burst 
of activity, which expels the accumulated energy, ensues. This process can be 
renewed, much like a domino effect, leading to a large transport event. Spe- 
cific conditions of the thermal convection in the container make these states 
metastable and produce random reversals. 

2 Multiscaling properties of the mean wind and SOC 

The probability density function (PDF) of the life-times of the metastable 
states, T, observed in the laboratory data on the wind, exhibit the character- 
istic power-law 

p(t) ~ r-\ (1) 

as shown in figure 1. 

Different physical mechanisms can lead to this power law. It could be the "ex- 
change of stabihty" under turbulent noise effect [12], or it could be the SOC 

as mentioned in the Introduction. These mechanisms differ substantially. The 
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Fig. 1. The probability density function of the hfe-times of the metastable states of 
the wind (in log-log scales). The dashed line indicates the power law (1). 

differences can be ascribed in part, almost trivially, to the lack of long time 
correlations for the avalanches in the "monoscale" SOC models, leading to 
the inability of SOC to incorporate the far-from-equilibrium characteristics of 
hydrodynamic systems such as turbulent convection. This trivial distinction 
between the two mechanisms disappears if we consider the two-dimensional 
Bak-Tang-Wiesenfeld (BTW) model [18], [19], which obeys a specific form of 
multiscahng for the PDF of several avalanche measures [20], [22], [23]. This 
model has bursts with long-time correlation, and other turbulence-like in- 
termittent properties, if studied on the time scales of its waves (see below). 
The question now is this: Can one still distinguish between turbulence and the 
multiscahng SOC in this situation? Or, are there still significant differences be- 
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tween the model and the specific physical phenomenon? We shall address this 
issue briefly in this paper, using as example the life-time r of the metastable 
states of the mean wind. 

The BTW prototype model of self-organized criticality is defined on a square 
lattice [19]. The number of "grains" stacked on a given site i is denoted by 
Zi — 0, 1,2, . . . If Zi < 4, Wi, the configuration is stable by construction. A 
random site k is selected and a grain is added to it, thus increasing Zk to 
Zk + 1. If, in the process, Zk > 4, immediately Zk is reduced to Zk — 4. This is 
called "toppling" . The expelled four grains are received one each by the four 
nearest neighbor site of k. The grains disappear if /c is a boundary site. It is 
clear that the toppling at site k may lead to topplings at the next microscopic 
time step, and so on. Thus, an avalanche made by a total number s > 
of topplings occurs before a new stable configuration is reached, and a new 
grain is added to a random site. After many additions the system reaches a 
stationary critical state in which the properties of the avalanches are sampled. 

A key notion of the BTW model is that of wave decomposition of avalanches 
[26] . The first wave is obtained as the set of all topplings which can take place 
as long as the site of addition is prevented from a possible second toppling. 
The second wave is constituted by the topplings occurring after the second 
toppling of the additional site takes place and before a third one is allowed, 
and so on. The total number of topplings in an avalanche is the sum of those 
of all its waves. The wave size, s, has the power-law PDF 

Pis) ~ s-\ (2) 
Returning now to convection, as mentioned in the Introduction, the times 



5 



42 r 




In t 

Fig. 2. Moments of the "local dissipation rate" of the metastable states duration 
times (or, in this case, their sizes). The straight lines (the best fit) are drawn in 
order to indicate (in log- log scales) the multiscaling given by (4). 



taken to switch from one metastable state to the other are very short in com- 
parison with the interval between switchings; in addition, the magnitude of the 
mean thermal wind velocity is approximately constant during the metastable 
events. Therefore, the size of the events in the experiment (in SOC terms) is 
proportional to their duration; that is, s ~ r (compare equations (1) and (2), 
[27]). 

With this identification, one can take the analogy a step further. A character- 
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Fig. 3. Multiscaling exponents from (4) for the mean wind (circles) and for the 
two-dimensional BTW prototype model of SOC (triangles). 



istic similar to local dissipation rate in turbulence, 

t 

et^J2i^k+i-Skf/t, (3) 

k=l 

was introduced in the Ref. [23] in order to characterize multiscaling properties 
of SOC. The multiscahng (if it exists) has the form 

Such multiscaling was observed in Ref. [23] for the two-dimensional BTW 
prototype model of the SOC. Since s ~ r in our case (see above and [27]), we 
calculated the analogous "local dissipation rate" for the duration r obtained 
in our experiment. The multiscaling and corresponding scaling exponents /Xp, 
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given by (4), are shown in figures 2 and 3, respectively. Circles in figure 3 
correspond to the exponents computed for our data whereas the triangles 
correspond to computations performed in [23] for the two-dimensional BTW 
model of SOC. 

3 Concluding remcirks 

The fluid flow here has many features such as jets, boundary layers (both 
attached and detached), plumes, and so forth. These are hydro dynamical en- 
tities which, in the end, may have to be explained through their hydrodynamic 
origin. Amidst this complexity, however, the comparison of the multiscaling 
exponents shown in figure 3 suggests that the multiscale SOC is a possible 
model for the observed reversal of the thermal mean wind. We might ask: Is 
the mechanism in the fiow indeed the same as that of two-dimensional BTW, 
or does the agreement of multiscaling exponents /ip shown in figure 3 char- 
acterize some class of the multiscale SOC universality? In some sense, the 
first question looks for more detailed dynamical understanding, while the lat- 
ter looks for a statistical analogy. We cannot answer the first question with 
any certainty yet, but the answer to the second question seems to be in the 
affirmative. 
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